Spin current injection by intersubband transitions in quantum wells. 
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We show that a pure spin current can be injected in quantum wells by the absorption of linearly 
polarized infrared radiation, leading to transitions between subbands. The magnitude and the 
direction of the spin current depend on the Dresselhaus and Rashba spin-orbit coupling constants 
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O , and light frequency and, therefore, can be manipulated by changing the light frequency and/or 

applying an external bias across the quantum well. The injected spin current should be observable 
either as a voltage generated via the anomalous spin-Hall effect, or by spatially resolved pump-probe 
optical spectroscopy. 
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Spin current is an interesting physical phenomenon in its own right, and could have application in the delivery 
and transfer of electron spins in spintronics devices. From a fundamental point of view, various issues raised in the 
theory of this effect are far from being satisfactorily settled. As was shown by Rashbai, a spin current exists even 
in the equilibrium state of a two-dimensional (2D) electron gas with spin-orbit (SO) coupling. The application of an 
external electric field has been suggested as a strategy for driving the system out of equilibrium and inducing a spin 
current exhibiting transport effects. Mal'shukov et al£ and Governale et al& suggested applying a time-dependent 
bias across a semiconducting heterostructure, thus modulating the strength of the SO coupling and generating a spin 
current. Murakami et alfi and Sinova et al~ have shown that an in-plane electric field can cause a spin current, 
leading to the "intrinsic spin-Hall effect" . Another possibility for the injection of spin current is coherently controlled 
' optical excitations between the valence and the conduction band, as predicted by Bhat and Sipe§*£ and observed 
■ experimentally in bulk crystals^ and quantum wells (Q Ws)AS. 

Here we show that a spin current can be injected in QWs by infrared (IR) light absorption, driving transitions 
between different subbands. The injection of spin-polarized electric current in QWs due to intersubband transitions 
caused by circularly polarized radiation has already been observed by Ganichev et al^ 1 . In contrast, here we investigate 
' a pure spin current, where electrons moving in opposite directions have opposite orientations of spins, not accompanied 
by a net electrical current. We show that the strength and direction of this pure spin current can be manipulated by 
modulating the SO coupling strength via applied bias 12 and/or adjusting the light frequency. 

As an example we consider the (Oil) GaAs QW, where the electron spins have a considerable out-of plane com- 
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ponent, thus making possible the observation of the pure spin current by detecting the voltage generated via the 



anomalous spin-Hall effec1ii2i±i. The first two subbands in the well are typically separated by the energy Ku>o w 100 
meV; the exact value depends on the width of the QW, dopant concentration, and the boundary conditions. The SO 
J-^ ' Hamiltonian for the (Oil) QW, Hso — Hjj + Hr, is the sum of a Dresselhaus termi^, Hd, originating from the unit 
cell inversion asymmetry, and a Rashba termi&, Hr, originating from the asymmetric doping and/or a bias applied 
across the well: 



H% ] = a [ Z ] k y a°F H (k), H% ] = (a*k y - a^k x ) , (1) 

where n is the subband index, k = (k Xl k y ) is the in-plane wavevector of the electron envelope function, Fi n ] (k) =1 — 
(ky — 2fc 2 ) A 2 n j , where A[„] depends on the QW width w, and the a 1 are the Pauli matrices. The z— axis is perpendicular 

to the QW plane and the in-plane axes are: x — [100] and y — [011]. The parameters ajp and depend on n; in 
the model of rigid QW walls one has = —a n 2 {-k/w) 2 /2, where «o is the Dresselhaus constant for the bulk, and 
Ar„i = w/nir 15 . The deviation of .FLi(k) from unity becomes important at electron concentrations N c \ ~ 10 12 cm" 2 . 
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The spin-related energy is given by -Ego (k) = ^ (a^' fc a i<]„](k)^ + ^a^'fc^ , with " up" (u) and " down" (d) states 
having energies (k) = ±-Ego (k), and leads to the subband spectra: 



^±4 (k), Es3 = ^ + ^±E s 2 o (k). (2) 



where m is the electron effective mass and the indices Si, S2 describe the u(+) and ef(— ) spin states in the subbands 
n = 1 and n — 2, respectively. The corresponding spin eigenstates result in expectation values of the spin 
components: 

[n] , p [n] 

(KW Z \K) = ± D I*™, > 1^") = ±^n7r(^-M, (3) 

^so W -^so W 

where upper(lower) sign corresponds to the u(d) state and a« = (a x ,a y ). 

There is not yet consensus in the literature on the fundamental description of spin current, and the effect of disorder 
on it, as discussed e.g., in Reftf; spin current is not a "true" current, in that its density does not satisfy a continuity 
equation describing the evolution of a spin density^. Nonetheless, we introduce a "physical" definition of spin current 
per electron as: 

j?(k,s n ) = J-<eK^ + <A*K"> ; (4) 

where fi and (3 are Cartesian indices. Velocity components v% = dH/frdki are the sums of normal i^n = hki/m and 
anomalous terms given in our model (Eq.(l)) by: 

[n] [n] [n] 

4tn = -^-o* + 4*^^o*Af B]l 4"L = ^ (l - (3*J - 2*2) Af n] ) a*. (5) 

Below we consider only the spin current components associated with the z— axis spin projection. First we calculate 
the equilibrium spin current at typical experimental conditions, where only the first subband is occupied, and then 
find the changes induced by the intersubband excitations. For this purpose we introduce the equilibrium Fermi 
distribution function for two spin projections in the first subband: 

/±( k ) = FTn \ 1 > (6) 



exp 



(h 2 k 2 /2m± E l sl(k) - fj\ /k B T 



where fi is the chemical potential for a given N e i, and ks is the Boltzmann's constant. The spin current density 
component J* is the sum of the normal, J* n , and the anomalous, Jy an , parts. By integrating jj- (k,si) over the 
equilibrium state we obtain: 

d 2 k I, /• .... , </-/,• ] 




(/+ (k) + /_ (k)) (1 - (3fc^ - 2k 2 x ) A?) —2 +-/(/+ (k) - /_ (k)) fc 



(7) 

where (</>k 1 |°' 2 \4>k ) 1S defined in Eq.(3), and J* — by symmetry. The contributions Jy an and J* (first and 
second term in Eq.(7), respectively) almost cancel each other. At T = each of them is close in absolute value 

to iV e i/2 ^1 — 7riVeiA^ l ]/2Y and Jy/Jy <an ~ (ma$ /tfkpj <C 1 where Tikp is the Fermi momentum (see also 
Rashbai). We show Jy/Jy aa as a function of N e i in Fig. 1. 

Now we can investigate the spin current injection by linearly-polarized IR radiation due to the intersubband tran- 
sitions, as shown in Fig. 2a. The external field is a pulse E(t) = £ (t) exp (— iu>t) + c.c. with the carrier frequency 
lu, and slowly varying amplitude £(t) of duration r. We consider oblique incidence with £(t) lying in the plane of 
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incidence. The radiation frequency u> is close to uiq, with a detuning Q = uj — ujq, such that it can cause transitions 
between the subbands, with TiSl being of the order of few meV. For r ^> uj^ 1 the exact shape of the pulse has no 
influence on our results; however, to have the possibility of momentum-selective excitations as shown in Fig. 2a one 
needs sufficiently long pulses, with r > Ti/aokF ~ 1 ps, for old ~ 10~ 9 eV-cm, a typical value of the Dresselhaus 
coupling 1 ™. This condition also implies applicability of Fermi's Golden Rule, since the pulse contains many periods of 
the field oscillations. Since , and, in turn, the spin states and anomalous velocities depend on the subband, 
the intersubband transitions can cause the injection of a spin current. The ratio a^J p r n ](k)/a^, which determines 
the direction of the effective SO field acting on the spin, depends on the subband. Therefore, the spin states in 
different subbands are not mutually orthogonal, so (cj)^ 1 </>£. 2 ) ^ 0, and, "spin-flip" transitions u < — ► d are allowed 
with linearly polarized IR light absorption. The transitions <f>^ — > (f>^ occur in the vicinity of the resonance curves 
in the momentum space, determined by the k = k£ 2,Sl (r2) where k* 2,Sl (ri) is specified by the constraint of energy 
conservation. For a given Q there are in fact two such curves. In our case E§> (k) > E 1 ^ (k) for all k. Therefore, for 

> the transitions d — > u and u — > u are allowed, while for Q < we obtain d — ► d and u — > d transitions. As one 
can see in Figs. 2a and 2b, fc* 2,Sl (0) is larger for the "spin-conserving" than for the "spin-flip" -transitions. 

The transition matrix elements depend on the spin states in both subbands, and can be factorized in the dipole 
approximation as: 



M(k s /' sl (n))=Se 




where 6o an d 0i are, respectively, the incidence and refraction angles, sin#i = smOo/y/e, e is the dielectric constant, 
e is the electron charge, and (p^'(z), ip^(z) are the envelope electron wavefunctions in the subbands n — 1 and n = 2, 
respectively. A transfer of one electron to the second subband injects a spin current: 

Af y (k;s 2 , Sl ) - f y (k, S2 ) - j* (k, Sl ) , (9) 

where we neglect the small photon momentum. The incident radiation injects the concentration of electrons in the 
second subband JVja, with a rate dN% (f2) jdt and, correspondingly, drives the spin current density component with the 
rate dAJ* (f2) jdt. The injection rates can be written as: 



dAJ*(») _ h , dN a (0) _m {S) (10) 

dt 2 dt ' dt Tito 

where £(f2) characterizes the effective speed of electrons forming the pure spin current, (S) = (c/2tt)£ 2 is the radiation 
power per unit area, and £(f2) is a dimensionless function. 

Within Fermi's Golden Rule the speed characterizing the spin injection is obtained as: 



E / /a t (k) I I O I 2 Aj* (k; S2 , ai ) dk s r ^ (n)/v s k ^ 

k(n) = ^ ? , (ii) 

E //,x(k)KCIOI 2 *r 2 ' SI (W sl 

si ,s 2 J 

with the velocity associated with the joint density of states given by: 

< 2 ' sl = ?k^-^- (12) 

The integration in Eq.(ll) is performed along the resonance curves. With the increase of |f2|, kf. 2 ' Sl (^l) increases and 
eventually arrives at regions of small electron occupancy, as can be seen in Fig. 2b. Hence, dAJ* (fi) /dt and dN2/dt 
become small at 7i|fi| larger than some critical h£l c (a few meV) determined by the condition min fc* 2 ' Sl (f2 c ) > ko, 
where kg — kp or fcp = y/mksT/h in the degenerate and non-degenerate gas, respectively. 
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The photoinduced spin current is the sum of normal A J* n and anomalous A J* an contributions, each containing 
spin-flip (si S2) and spin conserving (si = s 2 ) terms. The anomalous spin-conserving term is of the order of 
(a^} — ajp^ N 2 , while the other terms depend on the difference of the ratio ~~ a ^R / a ^D an< ^ \n]- An estimate 

of the relative contributions is: 

A J* n (si = 82) _ fiV 
AJ^ an ( s i = s 2) ma D 

Due to a large prefactor 7i 2 kp /map, which is the ratio of the normal and anomalous velocities, the normal term can 
be large and lead to a change in the sign of the spin current at particular light frequencies, as seen in Figs. 3a and 
3b. In Fig. 3a we present the speed C(f2), while in Fig. 3b we show the normal and anomalous parts of the injected 
spin current density. The spin-flip contribution in both the normal and anomalous terms is much smaller than the 
"spin-conserving" one. Recently, GolubAS demonstrated that the direction of electric current induced by interband 
light absorption in QWs can depend on the light frequency. In his scenario the change occurs as new subbands are 
accessed, and thus appears on a scale of 100 meV. In our scenario for pure spin current injection, the change occurs 
on a much smaller scale. 

Now we estimate the magnitude of the injected spin current assuming that the contributions of the anomalous 
and normal terms are of the same order of magnitude. Fig. 4 presents the efficiency of the energy absorption £(f2) 
(Eq.(10)). The concentration of the electrons excited to the second subband can be estimated from Eqs.(8)-(10) as 
N 2 ~ 2Tr(e 2 w 2 k F /e 2 hca D )(S)T. At e = 12, k F ps 10 6 cm" 1 ,!!) = 100 A, 9 close to tt/4 and N B \ ps 10 12 cm" 2 we 
obtain: A/2/iVei ~ 10~ 6 ^(S')/(W/cm 2 )^ • (r/ps). Under excitation of a 1% fraction of electrons, achieved at (S) ps 10 

kW/cm 2 and r ps 1 ps, the corresponding effective current density eAJy/h ps 1 mAmp/cm. This is of the same 
magnitude as would be generated by the ac spin pumping in the n = 1 subband, as proposed by Mal'shukov et al£, 
but the effect would operate on a nanosecond time scale, as opposed to the picosecond time scale relevant here. 

Having found the magnitude of the spin current, we discuss its experimental observation. A possible technique 
is the measurement of the voltage generated by the anomalous spin-Hall effect due to scattering of electrons by 
impurities. The spin current AJ* causes a spin-Hall bias V s u along the x axis. Its magnitude can be estimated as 
V s h ~ t&n(9 S H )V e s , where 9 s h is the spin-Hall angle and V c fi is the effective lateral bias that would cause a current 
density eAJy/h. As follows from the discussion preceding Eq.(13), the corresponding current density is of the order of 
eN 2 {ap>/h). The bias V c g that would cause this current density is: V g ps L {N 2 /N c \) au/KfJ,, where /i is the mobility, 
and L ps 1 cm is the lateral size of the system. At \i ps 10 5 cm 2 /(Vs), N2/N c \ ps 10~ 3 , and an/h ~ 10 6 cm/s we 
obtain: V e s/L ps 10 -2 V/cm. The spin-Hall angle was estimated by Huang et ali& as 9 s h ~ 10 -3 , which would 
lead to V s h ~ 10~ 5 V. Their model assumed charged dopants embedded directly in the QW, which considerably 
overestimates the magnitude of the effect when only a remote doping is present. For this reason, 10~ 5 V is clearly an 
upper estimate of the spin Hall bias. Nonetheless, even a bias smaller by two orders of magnitude than this would be 
experimentally accessible^*. 

Another possibility for observing the pure spin current is spatially resolved pump-probe spectroscopy, as applied 
by Hiibner et al. 9 and Stevens et alML to investigate the spin current injected by interband transitions. In those 
experiments the centers of the spin-up and spin-down of excited electron distribution were separated by approximately 
20 nm. In the experimental situation considered here, the spin-polarized spots can be separated by distances of the 
order of the electron free path I s» (fikp jra) Tfe, with Tk being the momentum relaxation time. At mobility /i ps 10 5 
cm 2 /(Vs), one obtains I ps 10 3 nm, and so a possible approach would be to observe this separation experimentally 
by using a linearly polarized IR light as a pump and circularly polarized light as a probe of the spin-dependent 
transmission. In a real sample, of course, we have to expect some inhomogeneity in the spin-orbit interaction due 
to quantum well thickness variations, dopant fluctuations, inhomogeneous strain, and the like^S. We are currently 
investigating the consequences of such inhomogeneity, and will return to it in a later communication. 

To conclude, we have shown that a pure spin current can be injected in QWs by IR intersubband absorption, 
calculated its magnitude, and found that it could be measured experimentally. The dependence of the spin current 
on the light frequency, and on the Rashba SO coupling parameter, opens the possibility of its manipulation applying 





(13) 
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an external bias and by changing the light frequency. The spin current should be observable by anomalous spin-Hall 
effect measurements or by pump-probe optical spectroscopy. 
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FIG. 1: Spin current Jy/Jy^ n as a function of the electron concentration N e \. Dashed curve : Ai = 0, solid curve Ai = to/ 
the QW width w = 80 A. The parameters are: «p =-0.3xlCP 9 eVcm, a^=-0.3a^ , ksT— 25 meV, m = 0.066mo, where mo 
a free electron mass. 
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FIG. 2: (a) The intersubband transitions leading to the injection of pure spin current. Thick arrow line corresponds to SI = 0. 
Thin arrow lines correspond to transitions at S7 > 0. (b) Resonance curves k = k° 2,sl (Sl). Solid (dash) lines describe the 
spin-conserving (spin-flip) transitions. In each case the outer curve is for TiQ. — 2 meV and the inner curve for TiQ. = 1 meV. 
The circle marked as &f is the Fermi line at N c \ — 10 cm" . The parameters are: a ™=-0.3xl(r 9 eVcm, a£ ] =-0.3a£ ] , 
a^=4a^, a^=-0.5a^, and Ar n i = w/nn. 
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FIG. 3: (a) The speed ({CI) f° r different Rashba coupling constants a R (values in 10~ 9 eVcm units are presented near the 
curves). Other parameters are the same as in Fig. 2, N e \ — 10 12 cm -2 , and ksT— 25 meV. The direction of the spin current 
can be altered by changing the Rashba parameter, (b) Components of the induced spin current as the function of the photon 
frequency for a 1 '^ = 0.6 x 10 -9 eVcm. 




FIG. 4: ^-dependence of for different Rashba parameters q^'. The units of incident light power density are W/cm 2 . 
take an incidence angle 9o = 7r/4, e = 12, and w = 80 A. The infinite barrier approximation is used for the calculation. 



